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^ ■ The canonical Monte-Carlo is used to study the phase transitions from the 

low-temperature ordered phase to the high-temperature disordered phase in 
the two-dimensional Falicov-Kimball model with correlated hopping. As the 
low-temperature ordered phase we consider the chessboard phase, the axial 
striped phase and the segregated phase. It is shown that all three phases 
persist also at finite temperatures (up to the critical temperature r c ) and that 
the phase transition at the critical point is of the first order for the chessboard 
and axial striped phase and of the second order for the segregated phase. In 
addition, it is found that the critical temperature is reduced with the increasing 
amplitude of correlated hopping t' in the chessboard phase and it is strongly 
r**~ ■ enhanced by t' in the axial striped and segregated phase. 
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1 Introduction 

The Falicov-Kimball model is one of the simplest yet most versatile models of strongly 
correlated electron systems on the lattice pQ. The model describes a two-band system 
of localized / electrons and itinerant d electrons with the short-ranged f-d Coulomb 
interaction U. The Hamiltonian is 

Ho = E tijdfdj + UJ2 tifidfdi + E f £ f+f t , (1) 

ij i i 

where f^~, f\ are the creation and annihilation operators for an electron in the local- 
ized state at lattice site % with binding energy Ef and df, di are the creation and 
annihilation operators for an electron in the conduction band. The conduction band 
is generated by the hopping matrix elements ty, which describe intersite transitions 
between the sites i and j. Usually it is assumed that t^ = — t if i and j are nearest 
neighbours and ty = otherwise (the conventional Falicov-Kimball model). 

The model has been used in the literature to study a great variety of many-body 
effects in metals, of which valence and metal-insulator transitions, charge-density 
waves and electronic ferroelectricity are the most common examples [21 El IH IS]- 
It has been applied to a variety of lattices, one [HI U\, two [SI [SI [TO], three [TTj . 
and infinite dimensional [12], and occasionally to small clusters [T3| O [TS]. Exact 
results are available in very few instances [T2| [TBI El 02] and general theorems 
have been proved for special cases [HI US]- In spite of the existence of an analytic 
solution in d = oo dimension [T2| I2U] and an impressive research activity in the 
past, the properties of this seemingly simple model are far from being understood, 
especially for nonzero temperatures. For example, it is well known [2TJ that the 
ground-state phase diagram of the model exhibits a rich spectrum of charge ordered 
phases including various types of axial and diagonal striped phases, the chessboard 
phase, the segregated phase, etc., but only a little is known about the temperature 
stability of these phases [221 E2]- Similarly, only a little is known about the type 
of phase transitions from the ground-state ordered phases to the high-temperature 
disordered phase [21]. From this point of view the most explored phase from the 
above mentioned ones is the chessboard phase that is the ground state of the model 
at the half-filled band case (Ef = 0, Nf = Nd = L/2, where L is the number of 



lattice sites). For this case there exists the exact proof [9 J of existence the phase 
transition from the low-temperature ordered phase (the chessboard phase) to the 
high-temperature disordered phase at finite critical temperature r c (for dimensions 
d > 2) that strongly depends on the local Coulomb interaction U. In addition, the 
numerical simulations within the grand-canonical Monte-Carlo showed that the phase 
transitions are of the first order for small and intermediate values of the Coulomb 
interaction U and of the second order for strong interactions [2H [25]. In the current 
paper we extend the numerical study of the temperature induced phase transitions 
also on the case of phase segregated and striped phases. Moreover, we consider here 
a more general situation 

H = H + H« (2) 

with the correlated hopping term 

H^t'Y^Uth + ftfMfd^ (3) 

(hi) 
that represents a much more realistic description of electron hopping in rare-earth 
compounds [26J. As was shown in our previous papers [27J [28] this term has strong 
effect on the formation of charge ordering in the ground state and therefore it should 
be taken into account in the correct description of electronic correlations in rare- 
earth materials. In particular, we have found that already relatively small values of 
the correlated hopping term lead to a stabilization of new types of charge ordering, 
even in the half-filled band case, where the ground state at if — is the chessboard 
phase for all nonzero U. The comprehensive ground-state phase diagram of the 
two-dimensional half-filled Falicov-Kimball model with the correlated hopping in the 
if — U plane is presented in [27]. It consists of three different phases, and namely 
(i) the chessboard phase located in the central region of the phase diagram along 
the U axis, (ii) the axial striped phase located below (U > 0) and above (U < 2) 
the chessboard phase, and (iii) the segregated phase located above the axial (U < 2) 
and chessboard (U > 2) phase. Since these phases represent the most prominent ex- 
amples of charge ordering observed experimentally in strongly correlated materials, 
like cuprates, nickelates and cobaltates, we have decided to perform exhaustive nu- 
merical studies of the half-filled Falicov-Kimball model with the correlated hopping 
with a goal to answer the questions about the temperature stability of these phases 



and the type of phase transitions from the low-temperature ordered phases to the 
high-temperature disordered one. 

2 Method 

Since in this spinless version of the Falicov-Kimball model with correlated hopping 
the /-electron occupation number ff fi of each site % commutes with the Hamiltonian 
(2), the /-electron occupation number is a good number, taking only two values: 
w{ — 1 or 0, according to whether or not the site i is occupied by the localized / 
electron. Therefore the Hamiltonian (2) can be written as 

(4) 





H — y^ hij{w^)dfdj 




{i,3) 


where hij{w*) - 


= tij(wf) + Uw{5ij and 




tij{w f ) = Uj + lfij(w{ + W f j 



(5) 

Thus for a given /-electron configuration w^—{w{,W2,- ■ ■ ,w£}, defined on the two- 
dimensional lattice of L sites, the Hamiltonian (2) is the second-quantized version 
of the single-particle Hamiltonian h(w'), so the investigation of the model (2) is 
reduced to the investigation of the spectrum of h for different configurations of / 
electrons. Since we are interesting in the half-filled band case, where both the total 
number of / and d electrons are fixed to L/2, the numerical calculations at nonzero 
temperatures are done exclusively in the canonical ensemble. In this formalism the 
partition function and the internal energy corresponding to the model Hamiltonian 
(2) can be written as: 



Z = E ^ /T , E = Y,e t (wf)wt (6) 

(E) = £ Ee- E '\ (7) 

where r = ksT and the summation goes over all possible L\/Nf\(L — Nf)\ distribu- 
tions w* of / electrons on L lattice sites and L\/N^.(L — AT d )! distributions w d of d 
electrons on L single-particle energy levels £$ corresponding to h(w*). In the next 
step the summation over all / and d distributions is replaced by the Monte-Carlo 
summation with the statistical weight e~ E l T jZ. 



To identify the transition temperatures from the low-temperature ordered phases 
to the high-temperature disordered phase and the type of the phase transition we 
have calculated numerically the specific heat C = {(E 2 ) — (E) 2 ) / (Lt 2 ) , the thermal 
average of the /-electron occupation w s = (it)*) and the energy distribution P(E). 
The numerical calculations are done exclusively at U = 0.5, since the ground-state 
phase diagram exhibits the richer spectrum of solutions in the weak and intermediate 
coupling regions in comparison to the strong coupling limit. 

3 Results and discussion 

To verify the ability of our method to describe the phase transitions at finite temper- 
atures we have started with the conventional two-dimensional Falicov-Kimball model 
(f = 0) at half-filling. As was mentioned above, the physical picture of temperature- 
induced phase transitions within this relatively simple model is well understood at 
present. For all finite Coulomb interaction U > the ground state of the model 
is the chessboard phase that persists up to critical temperature t c (U), where the 
system undergoes the phase transition to the homogeneous phase. The phase tran- 
sition is of the first order for U < 1 and of the second order for U > 1 [21]. Our 
numerical results obtained within the canonical Monte-Carlo method for C, w s and 
P(E) fully confirm this picture (see Fig.l). The specific heat curves exhibit a sharp 
low-temperature peak at r c ~ 0.028 that is connected obviously with the phase tran- 
sition from the chessboard phase to the homogeneous phase, as can be seen from the 
behaviour of the average /-electron occupation w s for temperatures slightly lower 
or slightly higher than r c . Moreover, the energy distribution function P(E) exhibits 
an apparent two-peak structure near the critical point r c (it can be considered as 
a superposition of two Gaussians), what in accordance with the theory of Challa, 
Landau and Binder [29J points on the first order phase transition at r c . 

Let us now discuss how this picture is changed when the correlated hopping term 
is added. Firstly, we have examined the case of small values of \t'\ for which the 
ground state of the model is still the chessboard phase [2Z] • The typical examples of 
C, w s and P(E) from the positive and negative region of t' are displayed in Fig. 2 
and Fig. 3 for t' = —0.3 and t' = 0.3. One can see that the correlated hopping term 



(in the limit of small \t'\) does not change qualitatively the picture of temperature 
induced phase transitions found for t' = 0. For both, positive and negative if, there 
is the first order phase transition from the low-temperature ordered phase to the 
high-temperature disordered one, similarly as for t' = 0, and the only difference 
between these cases is that the correlated hopping term reduces slightly the critical 
temperature r c of the phase transition. 

Therefore, in the next step we have turned our attention to the physically much 
less explored type of configurations, and namely, the axial striped configurations that 
are ground states of the Falicov-Kimball model for the intermediate values oft' (|t'| ~ 
0.5). Note, that for the axial striped phase even the fundamental question concerning 
the temperature stability of this phase has been not answered till now. This is caused 
by the fact that it is very difficult to find this phase in the pure form. For example, in 
the conventional Falicov-Kimball model {if = 0) the axial striped phases are stable 
for a relatively wide range of model parameters [21], but only in mixtures with 
other phases (e.g., the empty configuration). In addition, strong finite-size effects 
have been observed on the stability of these mixtures and therefore it is practically 
impossible to do any conclusions concerning their stability at finite temperatures 
from the numerical calculations on finite clusters. However, in the Falicov-Kimball 
model with correlated hopping the axial striped phase exists in the pure form for 
wide range of model parameters if and U, the finite-size effects on the stability of 
this phase at r = are negligible, and so the corresponding numerical study of the 
temperature stability of the axial striped phase can be performed straightforwardly. 

In Fig. 4 and Fig. 5 we present our canonical Monte-Carlo results for C, w s and 
P{E) obtained for two different values of if {if = 0.5 and t' = 0.55) from the region 
where the ground-state of the model is just the axial striped phase. Again, the specific 
heat curves exhibit the sharp low-temperature peak, the existence of which indicates 
the phase transition form the axial striped phase to the homogeneous phase. This 
was verified independently by calculating the average /-electron occupation w s and 
the energy distribution P{E) near the transition point r c , that clearly demonstrate 
the presence of the first order phase transition at r c . Since the critical temperature 
r c of the phase transition for both values of t' shifts to smaller values with increasing 
L, we have performed a detailed finite-size scaling analysis of the t c (L) dependence 
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to exclude a possibility that r c vanishes in the thermodynamic limit L — > oo. The 
resultant r c (L) dependencies are plotted as insets in Fig. 4 and Fig. 5. It is seen 
obviously that the critical temperatures r c for both t' = 0.5 and t' = 0.55 persist also 
in the thermodynamic limit, what means that the axial striped phase remains stable 
also at finite temperatures. In addition, our numerical results show that the critical 
temperatures for the axial striped phase are considerably higher in comparison to 
the critical temperatures for the chessboard phase. The same behaviour we have 
observed also for negative values of t' (£' = —0.7), however the critical temperature 
in this case was only slightly larger than one corresponding to t' = 0. 

With increasing if the half-filled Falicov-Kimball model with correlated hopping 
exhibits (at r = 0) the phase transition from the axial striped phase to the segregated 
phase [27] that takes place at t' ~ 0.6. Since the chessboard phase as well as the 
axial striped phase are both insulating and the segregated phase is metallic |27J, one 
can expect a fully different thermodynamic behaviour of the model for the last case. 
To verify this conjecture we have performed an exhaustive numerical studies of the 
temperature dependence of C, w s and P(E) for t' = 1. This study is important 
also from this point of view that the thermodynamic of the metallic phase has been 
examined till now only in a few cases [231 EH], while for the insulating phase (usually 
the chessboard phase) there is a number of analytical and numerical results [9| |3T| 132] . 

The results of our numerical calculations obtained for the specific heat C are 
shown in Fig. 6. To reveal the finite-size effects the calculations for C have been 
done on several different clusters of L=Q x 6, 8 x 8, 10 x 10, 12 x 12 and 16 x 16 
sites. We have found that the specific heat curves, in the low-temperature region, 
strongly depends on the cluster sizes, and therefore a very careful analysis has to 
be performed to find the correct behaviour of the model in the thermodynamic 
limit L — > oo. On small finite clusters (1 = 6x6 and L = 8 x 8) the specific 
heat exhibits only one-peak structure in the low-temperature region (r ~ 0.15). 
With the increasing cluster size L an additional peak is stabilized at slightly higher 
temperatures (r ~ 0.23), while the first peak is gradually suppressed and probably 
fully disappears in the thermodynamic limit. The behaviour of the average /-electron 
occupation shows (see Fig. 6) that the second peak in the specific heat corresponds 
to the phase transition from the low-temperature ordered (segregated) phase to the 



high-temperature disordered phase. 

The nature of this phase transition is, however, different in comparison to pre- 
vious cases. While the energy distribution function P(E) is double peaked for the 
chessboard and the axial striped phase near the transition temperature r c (the first 
order phase transition), P(E) exhibits the single-peak structure for the segregated 
phase, what points on the second order phase transition at r c . Comparing the ther- 
modynamic behaviour of the model in the chessboard, axial striped and segregated 
region one can find two other important differences, and namely, (i) the critical tem- 
perature of the second order phase transition is approximately ten times higher than 
the critical temperatures of the first order phase transitions, and (ii) the specific 
heat (in the low-temperature region) decreases exponentially for the chessboard and 
axial striped phase, while in the segregated phase the specific heat C{r) seems to 
have the linear behaviour indicating the Fermi-liquid behaviour for r < 0.08 (see the 
inset in Fig. 6a). The observation of the linear contribution to the specific heat in 
the low-temperature region (r < 0.08) is consistent with behaviour of the average 
/-electron occupation in this region (see Fig. 6c). One can see, that despite the in- 
creasing temperature (from to 0.08) the /-electrons preferably occupy only one half 
of the lattice leaving another part empty. Due to the on-site Coulomb interaction 
between the / and d electrons, the itinerant d electrons occupy preferably the empty 
part of lattice, where they can move as free particles yielding the linear contribution 
to the specific heat. 

In summary, we have studied the phase transitions from the low-temperature 
ordered phase to the high-temperature disordered phase in the two-dimensional 
Falicov-Kimball model with correlated hopping using the canonical Monte-Carlo. 
As representative examples of low-temperature ordered phases we have chosen the 
chessboard phase, the axial striped phase and the segregated phase. It was shown 
that all three phases persist up to critical temperature r c and that the phase tran- 
sition at the critical point is of the first order for the chessboard and axial striped 
phase and of the second order for the segregated phase. In addition, we have found 
that the critical temperature is reduced with the increasing amplitude of correlated 
hopping t' in the chessboard phase and it is strongly enhanced by t' in the axial 
striped and segregated phase. 



Acknowledgments 

This work was supported by Slovak Grant Agency VEGA under Grant No. 2/0175/10, 
Slovak Research and Development Agency (APVV) under Grant VVCE-0058-07. 
H.C. acknowledges support of Stefan Schwartz Foundation. H.C. thanks M. Zonda 
for stimulating discussion on the canonical Monte-Carlo. 

References 

[1] L.M. Falicov and J.C. Kimball, Phys. Rev. Lett. 22, 997 (1969). 

[2] D.L. Khomskii, Quantum Theory of Solids, edited by I.M. Lifshitz (Mir, Moscow 
1982). 

[3] T. Portengen, T. Ostreich, L.J. Sham, Phys. Rev. Lett. 76, 3384 (1996). 

[4] CD. Batista, Phys. Rev. Lett. 89, 166403 (2002). 

[5] P. Farkasovsky, Phys. Rev. B 77, 155130 (2008). 

[6] J.K. Freericks and L.M. Falicov, Phys. Rev. B. 41, 2163 (1990). 

[7] Ch. Gruber, D. Ueltschi and J. Jedrzejewski, J. Stat. Phys. 76, 125 (1994). 

[8] U. Brandt and R. Schmidt, Z. Phys. B 63, 45 (1986); 67, 43 (1987). 

[9] T. Kennedy and E.H. Lieb, Physica 138A, 320 (1986); E.H. Lieb, ibid. 140A, 
240 (1986). 

[10] Ch. Gruber, J. Iwanski, J. Jedrzejewski and P. Lemberger, Phys. Rev. B 41, 
2198 (1994); Ch. Gruber, J. Jedrzejewski and P. Lemberger, J. Stat. Phys. 68, 
913 (1992); T. Kennedy, Rev. Math. Phys. 6, 901 (1994). 

[11] R. Ramirez, L.M. Falicov and J.C. Kimball, Phys. Rev. B. 2, 3383 (1970). 

[12] U. Brandt and C. Mielsch, Z. Phys. B 75, 365 (1989); Z. Phys. B 79, 295 (1990); 
Z. Phys. B 82, 37 (1991); see also V. Janis, Z. Phys. B 83, 227 (1991); P.G.J. 
van Dongen and D. Vollhardt, Phys. Rev. Lett. 65, 1663 (1992); P.G.J, van 



Dongen, Phys. Rev. B 45, 2267 (1992); J.K. Freericks, Phys. Rev. B 47, 9263 
(1993); 48, 14797 (1993). 

[13] P. Farkasovsky and I. Batko, J. Phys.: Condens. Matter 5, 7131 (1992). 

[14] P. Farkasovsky, Phys. Rev. B 51, 1507 (1995). 

[15] P. Farkasovsky, Phys. Rev. B 52, R5463 (1995). 

[16] R. Lyzwa, Physica A 192, 231 (1993). 

[17] P. Lemberger, J. Phys. A 25, 715 (1992). 

[18] Ch. Gruber, J.L. Lebowitz and N. Macris, Europhys. Lett. A 21, 389 (1993); 
Phys. Rev. B 48 4312 (1993). 

[19] J.K. Freericks, E.H. Lieb and D. Ueltschi, Phys. Rev. Lett. 88, 106401 (2002). 

[20] J.K. Freericks and V. Zlatic, Rev. Mod. Phys. 75, 1333 (2003). 

[21] R. Lemanski, J.K. Freericks and G. Bannach, J. Stat. Phys. 116, 699 (2004); 
Phy. Rev. Lett. 89, 196403 (2002). 

[22] J. Wojtkeiwicz, G. Musial and L. Debski, phys. stat. sol. c 3, 199 (2006). 

[23] K. Czajka, M.M. Maska, phys. stat. sol b 244, 2427 (2007). 

[24] M.M. Maska, K. Czajka, Phys. Rev. B 74, 035109 (2006). 

[25] M. Zonda, P. Farkasovsky and H. Cencarikova, Solid State Communications 
149, 45 (2009). 

[26] J. Wojtkiewicz, R. Lemanski, Phys. Rev. B 64, 233103 (2001); Acta Physica 
Polonica B 32, 3467 (2001). 

[27] P. Farkasovsky and N. Hudakova, J. Phys.: Condens. Matter 14, 499 (2002). 

[28] H. Cencarikova and P. Farkasovsky, phys. stat. sol. b 242, 2061 (2005). 

[29] M.S. Challa, D.P Landau and K. Binder, Phys. Rev. B 34, 1841 (1986). 



10 



[30] H. Cencarikova and P. Farkasovsky, Czech. J. Phys. 54, D423 (2004). 

[31] P. Farkasovsky, Phys. Rev. B 54, 11261 (1996); Z. Phys. B 102, 91 (1997). 

[32] C.A. Macedo, L.G. Azevedo and A.M.C. de Souza, Phys. Rev. B 18, 184441 
(2001). 



11 



(a):f=0, U=0.5 



(b):f=0, U=0.5, L=16x16 




LU 



0.08 



0.06 
0.05 
0.04 
0.03 
0.02 
0.01 




x=0.028 
x=0.0285 
x=0.029 - 




(c):f=0, x=0.027 



w 
s 

• - 1 

• - 0.9 

• - 0.8 

• - 0.7 

• - 0.6 

• - 0.5 

• - 0.4 

• - 0.3 

• - 0.2 
■ - 0.1 

- 



(d):f=0, x=0.029 



Figure 1: The specific heat (a), the energy distribution (b) and the thermal average of 
the /-electron occupation (c-d) for the conventional Falicov-Kimball model (£' = 0) 
in two dimensions. 
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Figure 2: The specific heat (a), the energy distribution (b) and the thermal average 
of the /-electron occupation (c-d) for the two-dimensional Falicov-Kimball model 
with correlated hopping i! = —0.3. 
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Figure 3: The specific heat (a), the energy distribution (b) and the thermal average 
of the /-electron occupation (c-d) for the two-dimensional Falicov-Kimball model 
with correlated hopping i! = 0.3. 
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Figure 4: The specific heat (a), the energy distribution (b) and the thermal average 
of the /-electron occupation (c-d) for the two-dimensional Falicov-Kimball model 
with correlated hopping t' = 0.5. The inset shows the critical temperature r c as a 
function of the cluster size L. 
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Figure 5: The specific heat (a), the energy distribution (b) and the thermal average 
of the /-electron occupation (c-d) for the two-dimensional Falicov-Kimball model 
with correlated hopping t' = 0.55. The inset shows the critical temperature r c as a 
function of the cluster size L. 
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Figure 6: The specific heat (a), the energy distribution (b) and the thermal average of 
the /-electron occupation (c-f ) for the two-dimensional Falicov-Kimball model with 
correlated hopping t' — 1 . The inset shows the specific heat C in the low-temperature 
region for L = 16 x 16. 
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